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Abstract 

The dynamics of finite nonperiodic Toda lattice is an isospectral deformation of the 
finite three-diagonal Jacobi matrix. It is known since the work of Stieltjes that such 
matrices are in one-to-one correspondence with their Weyl functions. These are ra- 
tional functions mapping the upper half-plane into itself. We consider representations 
of the Weyl functions as a quotient of two polynomials and exponential representa- 
tion. We establish a connection between these representations and recently developed 
algebraic-geometrical approach to the inverse problem for Jacobi matrix. The space of 
rational functions has natural Poisson structure discovered by Atiyah and Hitchin. We 
show that an invariance of the AH structure under linear-fractional transformations 
leads to two systems of canonical coordinates and two families of commuting Hamilto- 
nians. We establish a relation of one of these systems with Jacobi elliptic coordinates. 
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1 Introduction. 



The Toda lattice is a mechanical system of iV-particles connected by elastic strings. The 
Hamiltonian of the system is 



N-l 



N-2 



Introducing the classical Poisson bracket 



N-l 



{f,g} = J2 



k=0 



d£dg_ _ df dg 
dq k dp k dp k dq k ' 



we write the equations of motion as 

Qk = ilk, H} = p k , 

p' k = { Pki H} = - e «fc-9*+i + e ,k - 1-,k , 



k = 1, 



N-l. 



We put g_i = — oo, = oo in all formulae. Following |TD|, [19] introduce the new variables 

D (/ fc -g fe+ l/2 



Cfc = e 1 



Vk = ~Pk- 



In these variables 



and 



N-l 



N-2 



A;=0 



k=0 



{Cfc, V k } = -Cfc/2, 

The equations of motion take the form 



{c k ,v k+1 } = c fe /2. 



'1.1: 



{v k , H} = c\ 



c* = {c fc , #} = c fc (^+i - v k )/2. 
These equations are compatibility conditions for the Lax equation L* = [A, L] , where 



v c 

C Vi Ci 








and 



2A 



■■ 

■■• 

c 

-c 

■■• 

••• 



CN-3 v N-2 CN-2 

Cat_ 2 V N -\ 



Ci 








-CN-3 C N -2 

-Cjv_ 2 
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The Lax formula implies that the spectrum Ao < • . • < Ajv-i of L remain fixed. It is known 

that the rational function w(X) = (R(X)8(0), 5(0)) where 



25 



since the work of Stiltjies, 
-R(A) = (L — XI)^ 1 is the resolvent, plays the key role in reconstruction of the matrix L from 
its spectral data. It was encountered later in the spectral theory of the Sturm-Liouwille 



operator, [p5[] , and received the name of Weyl function, 
continued fraction 



Simply by expanding it in the 



w(X) 



1 



A - v - 



A — v\ 



X-v 2 



r 2 
N- 



one can read the entries of L 
To da 



A - v N -i 

This fact is central to Moser's solution of the nonperiodic 



21 



At the same time Stiltjies method appears to be a computational trick and does 
not provide a conceptual explanation why the function w(X) actually determines L. 

An attempt to solve an inverse spectral problem for the finite Jacobi matrix using 
algebraic-geometrical approach was made at the beginning of eighties by McKean, [p0 |. It 
was realized that the corresponding spectral curve is a singular reducible Riemann surface. 

The recent interest in the spectral curves of finite Toda lattice stems from various sources. 
One is the work of Seiberg and Witten, [p23|, on sypersymmetric Yang-Mills theories. The 
smooth hyperelliptic spectral curves of the periodic Toda appear in the pure gauge N = 2 



SUSY Yang-Mills models in 4 dimensions |T2[ . The physically interesting limit of the theory 



which corresponds to transition from the smooth hyperelliptic curve to the singular one of 
2Cj was considered in ||. Another source of interest is the Camassa-Holm equation, [0]. Is 
is known, see G2|, 13] , that the dynamics of, so-called peakons solutions, is isomorphic to an 



isospectral flow on the space of finite Jacobi matrices. 

A solution of nonperiodic Toda within algebraic-geometrical approach was obtained re- 
cently in [pl|. The poles of w(X) determine the curve and the zeros specify the divisor of the 
Baker-Akhiezer function. This information uniqly specifies the BA function and therefore 
the matrix L. Whence that fact that w(X) determines the matrix can be considered as a 
consequence of the Riemann-Roch theorem which guarantees uniqness of the BA function. 
Using BA functions we obtained in the explicit formulae for the solution, the symplectic 
structure and two systems of canonical coordinates for it. 

The present paper combines ideas developed in 

of 



with algebraic-geometrical approach 



18 



It can be divided in two parts. In the first part of the paper (sections 2-5) we es- 
tablish relations between the standard objects of spectral theory and algebraic-geometrical 
constructions. We show how the BA function can be constructed from the orthogonal poly- 
nomial of the first kind and suitably normalized Weyl solution. The rational functions which 
map the upper half-plane into itself are defined by the formula 

Po Pn-i 



w(X) 



An -A 



+ 



Pk > o. 



1.2 



An-i — A 

They are parametrized by 2N parameters p's and A's. The rational functions corresponding 
to the Weyl functions of finite Jacobi matrices are specified by the condition Pk = 1- 
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We consider different reperentations of the Weyl functions. The first one, the ratio of 
two monic polynomials, 

leads to the standard Abel map. Another, the exponential representation, which employs 
Krein's spectral shift function £(z) 

MA) = -J_expE(A), 2(A) = f ^-dz 

A - A J z- X 

produces the Abel map in Baker's form. 

In the second part (sections 6-9) we study the Atiyah — Hitchin bracket on rational func- 
tions and its relation to the Hamiltonian formalism for the Toda lattice. As it was recently 
discovered by Faybusovich and Gehtman, || the AH bracket on rational functions can be 
written as 

{w(X),w(fi)} = y -^- (1.3) 

A — /i 

The bracket |TT] corresponds to the restriction of |1.3| on the 2N — 1 dimensional submanifold 
corresponding to Weyl functions of N x iV Jacobi matrices. We show that the restricted 
bracket has two systems of canonical coordinates on symplectic leafs of the foliation defined 
by level sets of the Casimir Afc- The existence of these two canonical coordinate sys- 
tems is a consequence of invariance of the AH bracket under the group of linear-fractional 
transformations 

, aw + b 



w — > w 



cw + d 



The first system of canonical coordinates is associated with iV poles of w(X) : Ao < 
. . . < Ajv-i- The half of variables is iV — 1 points of the spectrum A 1; . . . , Ajy-i- Another 
half are the functions of q(X) at these points 

o t = log tlh!M, * = i,...,*-i. 

<?(A ) 

This system is called action-angles coordinates. The verification of canonical relations is 
obtained by computing residues. The associated Hamiltonias are 



J 

They produce the standard Toda flows which preserve the spectrum. 

The second system is coming from N — 1 finite roots 71 < . . . < 'Jn-i of the equation 

= "T + • • • + T = 0. 

Ao - 7 Aat_i - 7 

The adjoint N — 1 variables are 

7r fc = log(-l)" + *p( 7fc ), k = l,...,N-l. 
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These are the divisor-quasimomentum coordinates. The Hamiltonians of this system are 

jY.~>- i 1 v '• 

These Hamiltonians produce the flows transversal to the isospectral manifolds. These flows 
preserve the divisor. 

Now we will explain the origins of the divisor-quasimomentum coordinate system. Jacobi, 
|13fl lecture 26, introduced his famous elliptic coordinates as finite roots 70 < 71 < . . . < 7jv-i 
of the equation 

W[V = 7 + • • • + T = 1. 

A - 7 Ajv_i - 7 

These N roots are considered to be the functions of iV independent p's while A's remain 
fixed. In the case of finite Jacobi matrices this choice of parameters 7 is bad. Due to the 
constraint Pk — 1 they are functionaly dependent. The same is true for any other constant 
instead of 1. Only for the special value the equation it; (7) = has one "unmovable" root 
70 = 00 and all other N — 1 finite roots are functionaly independent. Whence the divisor 
can be considered as a variant of the Jacobi elliptic coordinates. 

The transformation w(\) — > it/(A) = —l/w(X) defines the dual Weyl function it/. The 
roots of the equations 10(7) = become poles of the function it/ (A). The invariance of the AH 
structure allows to establish canonical character of the divisor-quasimomentum coordinates 
again by computing residues. 

We would like to conclude the introduction with the following remark. The tradional 
study of the integrable dynamics is a part of Hamiltonian mechanics with its' standard 
objects like Poisson and symplectic manifolds, vector fields, differential forms, etc. We 
demonstrate in this paper that the integrable dynamics on the space of Jacobi matrices or 
equivalently Weyl functions can be reformulated and studied purely in terms of complex 
analysis. 

Organization of the paper. In section 2 we review standard constructions of spectral the- 
ory of three-diagonal Jacobi matrices. In section 3 we introduce the spectral shift function 
and establish the trace formulae. In section 4 we consider the reducible Riemann surface and 
show how the Baker- Akhiezer function can be constructed using the orthogonal polynomial 
of the first kind and suitably normalized Weyl solution. In section 5 we consider representa- 
tion of the Weyl function as a ratio of two polynomials and exponential representation. We 
demonstrate how using these representations one can construct the Abel map in the ordinary 
form and in the Baker form, correspondingly. We also describe the range of the map and 
show its' one-to-one character. This section completes the description of the spectral theory 
and its' relation to algeraic geometry. In section 6 we introduce the Atiyah-Hitchin Poisson 
bracket and describe some of its' elementary properties. In the next, 7-th section we con- 
struct canonical variables for the AH bracket. We also construct the Dirac restriction of the 
AH bracket on the submanifod of Weyl functions. We introduce two systems of canonical 
coordinates for the restricted bracket in section 8. We also linearize the flow in terms of 
spectral shift function. The last 9 section describes the isospectral and transversal flows. 

Acknoledgements. The author thank H. Braden, H. McKean, S. Natanzon and I. Kriche- 
ver for stimulating discussions. 
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2 The Direct Spectral Problem. 



Most of the material of this section can be found in Jl| and presented here in order to set 
notations. Consider a finite Jacobi matrix 



v c 

C Vi Ci 








• • • Cat_ 3 V N -2 C N - 2 
• • • Cjv_2 V N -i 

It acts as a selfadjoint operator in the complex / 2 [0, N — 1] with standard orthonormal basis 
5{k) = (..., 1, .. .), k = 0, . . . , N — 1. The operator L has simple spectrum Ao < • • • < \n-i 

fc-th place 

corresponding to normalized eigenvectors e(Ajt) = (eo(Afc), . . . , ejv-i(Ajt)), = 0, . . . , N — 1. 
Let -E'(A) = Yl e (^fc) ® e (Afc) be an orthogonal spectral measure of L. Thus 

A fe <A 



XdE{z) 



and 



*w = (L-Azr l = /^ 

For any two vectors tt, u the Parseval identity holds 

5^(M,e(A fc ))(i;,e(Afc)) 



Mi" 



(2.1; 



(2.2) 



(2.3) 



We associate with L the eigenvalue problem 

voVo + c yt = Xy , 

coVo + viyi + cty 2 = Xyi, 

Cn-iy n -i + v n y n + c n y n+1 = Xy n , 

cn~2Vn-2 + vn-Wn-i + c N ^y N = \y N -x- 



n = 2 N-2: 



(2.4) 
(2.5) 

(2.6) 



The coefficient c^-i is defined by the formula c^-i = rifc=o c fc • For the system pT4| - pT6 
introduce the solution 



we 



P(A): 



and for the system [2.5| - [2T6 



Q(A): 



P_i(A) = 0, P (A) = 1, ...,P N {X); 



Qo(A)=0, Q 1 (X) = -,..., Q N (X). 
c 
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Let L[k, p ] be the truncated matrix 



Vk C k 
Ck V k+ i 






Cp-l 



Cp—1 "r, 



(2.7) 



Then for n = 1, 2, . . . ,N — 1 



p„(A) = (-ir det(L[0 LT 1] ~ A/) 



and 



g n (A) 



n-l 

n c fc 

k=0 



ir +l det ( L [l,n-l] ~ A/ ) 

' n— 1 

n 

A;=0 



(2.8) 



(2.9) 



Now we will use the solutions P and Q to give the formula [2.3| more concrete form. 



For P(A) = (P (A), . . . , Pv-i(A)) from we have 
e(A fc ) = P(A fe ) p k 



JV-l 

E Pn 2 (A.) 
n=0 



(2.10) 



Thus 2.3 takes the form 



where 



u, v) = I u(X)v (A) dcr(X) 



(2.11) 



u(X) = (it, P(A)), 5(A) = («, P(A)), da(A) = ^ 5(A - A fe )p*. 
Moreover, using |2.10| we have 



JV-l 



(e(\ k ),5(0)) = ^p- k >0 



(2.12) 



k=0 



This implies, in particular, that 5(0) is a cyclic vector for L. 
For the Weyl function defined as 

Qtv(A) 



Pv(A) 



formulae |2.8| - |2.9| imply 



w(X) 



-l^+MetCW-jj-AJ) 



(-l) 7V det(L - A/) 



AT-l 

■i)^ 1 n - a) 

s=l 

AT-l 

-irn(An-A) 

n=0 



(2.13) 
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where the roots A and 7 interlace 

A < 7l < -^1 < • • • < < 1n-i < Ajv_i 

due to the Sturm theorem. 



(2.14) 



By construction Qn + wPn = for all A. Formulae p^H2~6| produce (L — XI) (Q + wP) 
5(0) and Q + wP = R(X)5(0). The formula |J implies 



w(X) = (R(X)6(0),6(0)) 



(5(0),dE(z)S(0)) f da(z 



z-X 



z-X 



From 2.1 for the moments of the measure da we have 



s k = (L k S(0),8(0)) = / X k da(X). 



Using 



w 



(X) = -J2 s n X~ (n+1) , where s = 1. 



(2.15) 



(2.16) 



(2.17) 



n=0 



We conclude this section with derivation of trace formulae. To simplify notations we 



assume that A = 0. Then, 2.13 becomes 



w 



1 N ^ ( 

k=l v 



Ik - x 
x k -x 



(2.18) 



After simple algebra, 
Ik - A 



Xh — X 



00 \ 00 

*+E * xp~ -E A P^ where A ° = 1 - 

p=l p=0 



Fat her more, 



(A) = -E 



n=0 



AT-1 

PiH hpjv_i=n fc=l 



A" (n+1) . 



Comparing it with [2.17| we obtain the trace formulae 

AT-l 

«»= e ri A "- 



(2.19) 



PiH \-pN-i=n k=l 



The first few are listed below 

si 



k 

S2 = E A *+E A ^ A L> 



k\^k2 



S3 



E A *+ 2 E A ^ A U E A fci A fe2 A fc 3 ' etc. 

fc fci^fe ki^=k2^ks 
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To derive the standard trace formulae, from the resolvent expansion 

ILL 2 



R(X) = (L — XI) 



-1 



A A 2 A 3 



we obtain 



MA) 



(15(0), 5(0)) (L5(0),5(0)) (L 2 5(0),5(0)) 



A 

_I _ "Hi - v ° + c o 
A A 2 A 3 

Matching the coefficients in |2.17| and [2.2C| . 



A 2 



A 3 



(2.20) 



k 

E A *-E( A *) 2 > ete - 



3 The trace formulae via Krein spectral shift. 



We assume A = 0, then 2.13 becomes 



w(X) 



1 det(L[i iAf _i] - XI) 
Adet(L|KerL ± -A/) : 



(3.1; 



where L | Ker L 1 - ia a restriction of L on the orthogonal compliment to Ker L. By elementary 
transformations the ratio of two determinants can be put in the form 



w(X) 



-nf— 



--expS(A) 



A 



exp 



JV-l 

E 

s=l 



dz 
z- X 



exp / — - J : . 



A 



z — X 



(3.2) 



whereQ 

= n L\KerL^(z) ~ n L[1>N _ x] (z) 

is the Krein spectral shift function, [|l^] . This exponential representation of the Weyl function 
has much wider range of applicability then the formula |3.1| , which requires separate existence 
of determinants in the numerator and denominator. It can be obtained, for example, for 
infinite unbounded matrices under very mild condition on closeness of Lu and L \ Ker L^ . 



1 ul(z) is a counting function til(z) = ^{eigenvalues of L < z}. 
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One can obtain trace formulae in terms of = J z k ^(z)dz entering into the asymptotic 
expansion 

oo 

l+l) 

J n /x 

n=0 



h(a) = - x;/ n A-^ 



Expanding the exponent in 3.2 and matching the coefficients with 2.17 



si = -fo, 

c - ^ f 

52 - y-/i, 

P 

53 = fofi - f2 ~ etc. 

o 



Evidently, these formulae can be put in the form |2.19| using representation of as a 
difference of two counting functions. 



4 The Spectral Curve. The Baker- Akhiezer Function. 

The function w(X) determines the matrix L or in another words the functions w(X) = w(X, L) 
are coordinates on the space C of all N x N Jacobi matrices. The "index" A which labels 
the "coordinates" takes the values in C 1 \spec L. This statement goes back to Stieltjes, 



| 25|j . There are two standard ways to recover L from w. The first is to expand w(X) into 
continuous fraction, from which one can read the coefficients of L. The second is to construct 
polynomials orthogonal with respect to the spectral measure recovered from w(X). A three 
term recurrent relation for these polynomials is, in fact, the matrix L. Recently, the classical 
inversion problem, received a new, algebro-geometrical solution, fllSf. The main novel part 



of |18| is, so-called, Baker-Akhiezer function for a reducible curve. This construction is 
described below. 

We start with the standard Weyl solution Q + wP and note that ® +wP is a solution of 



2~5| - |2.6| which is equal to 1 at n = and vanishes at n = N for all A. The vector P is also a 
solution of |2.5| - |2.6| which is equal to 1 for n = and all values of A and vanishes at n = N 
for A = Afc. Thus we have a "gluing" condition: 

P(X k ) = ^±^(A fc ). (4.1) 
w 

In other words, at the points of the spectrum the function w(X) conjugates two solutions P 
and Q + wP which vanish at the left (n = — 1) or right (n = N) correspondingly. 

The singular algebraic curve T (Figure 1) is obtained by gluing at the points of the 
spectrum two copies of the complex plane. Define Baker-Akhiezer function on T by the 
formula, A = A(e): 

n 1 \ Q + wP/w(X) : if e G T_. 
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Figure 1: Riemann surface 



The function ip is continuous on T due to the gluing condition |4.1|. The BA function has 



the only simple poles at the points of the divisor (7^ — ), . . . , (7^-1, — )• At two infinities 
the BA function has poles of prescribed order. These data determine the BA function and 
therefore the operator L. Thus, the uniqueness statement can be viewed as a consequence 
of the Riemann-Roch theorem. This remark completes the description of the direct spectral 
problem. The inverse problem can be solved using an explicit formula for the time-dependent 
BA function (for details see [ITH] ) . 



5 The Abel map and the Jacobian. 

Let RatN be a set of all rational functions wich map the upper half-plane into itself, vanish 
at infinity and have N poles. Any function w(X) from Rat^ has the form 



N-1 



Afc — A 

fc=0 



with real poles at Ao < ... < Ajv_i and pk > 0. For A real below/above spectrum the 
function w(X) is positive/negative. Furthermore, 

and w(X) continuously changes from minus infinity to plus infinity, when A runs between 
two consecutive poles. Thus the function w(X) has exactly iV — 1 zeros 7's which interlace 
A's as in the formula [2. 14 . 
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Furthermore, any function from Rat^ can be represented as a ratio of two polynomials 

JV-l 

*(A) ^(-i)-n(T,-A) qoXN -. +qiXN -> 



w(X) 



p(X) 



N-l 

-1)^ n (An - A) 

n=0 



+ ■ ■ ■ + gjV-1 

X N + Po X N - 1 + ...+p N - 1 



(5.2) 



The polynomials are defined up to a multiple factor. In the formula [5.2| it is chosen such that 
the leading coefficient of the denominator is 1, similar to |2.13| . Evidently, the polynomial 
q(X) can be determined from its' values q(X ), . . . , g(Ajv_i) which are free parameters. 

Now we turn to the submanifold Rat' N with Pk — 1- These are the functions from 
Rat N which are Weyl functions of finite Jacobi matrices. The values q(X ), . . . , q(Xiy^i) are 
not independent anymore. Indeed, from the identity 



Pn 



p(X) ^ X n - X 
and condition p(X n ) = we obtain q(X k ) = p'(Xk)pk- Therefore, 

N-l 



~ X 



n=0 



p'(X r 



1. 



Due to the relation go = Y2Pk, an another way to say that w(X) G Rat' N is that the 
polynomial q(X) in [5[2] is a monic polynomial. 

For a function w(X) G Rat' N we define angle variables by the formula 



Ok = log • 



-l) h q(X k 
<?(A ) 



k = l N-l. 



(5.3) 



The are exactly k roots of q(X) between Ao and A& and it changes sign k times when A varies 
from Ao to A&. Whence variables #'s are always real. 

To clarify geometrical meaning of the variables #'s we introduce (Figure 2) a standard 
homology basis on the curve V corresponding to some Jacobi matrix. We define differentials 
ujk by the formula 



1 



z — Xk z — Ao_ 
Evidently, the normalization condition holds 



k 



N-l. 



u k = 2ni 5 p k , 



k,p = l,...,N 
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Figure 2: Basis of circles 



while 6-periods of cj's are real and infinite. For the variables 9, we have 



9 k = nik + log q(X k ) - log q(X ) 




■o 



dz 



(5.4) 



N-l y 

nik + } / 




Therefore, #'s are the values of the Abel map from the divisor 71, . . . < 7at-i into R N 1 , the 
noncompact real part of the Jacobian. It will be shown that this map is onto. 

Fix some matrix Lq. All matrices with the same spectrum Ao < ... < Ajv_i as Lq 
constitute a spectral class of L , which we denote by S(L ). The spectral class S(L ) is in 
1 : 1 correspondence with Weyl functions from Rat' N . 



Theorem 1 (i) The variables 7i, • • • , 7at-i are coordinates on S(L ). Any sequence of ' j's 
which occupies open segments Afe_i < 7fc < A^, k — 1, . . . , N — 1 corresponds to some 
matrix from S (L ) . 

(ii) The variables 9±, . . . , ^at_i are coordinates on S(L ). Any sequence of 9's from R N ~ 1 
corresponds to some matrix from S(L ). 



N-l 
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Proof, (i) The variables 7's determine the roots of q(X) and therefore the function w(X). 
Whence, 7's are coordinates. 

To prove that 7's are free pick any sequence of 7's and form q(X). Then, by Lagrange 
interpolation 



<?(A) 5p v( x n) 1 

r>l( 



PW t^P'^n)X n -X- 

It is easy to check all p n = q{X n )/p'(X n ) are strictly positive. It remains to prove that 
^2iPn = 1- Indeed, the formula 

^-2^ A x n) x-x n 

implies 

,. <?( A n) M A ) ?( A n) 



A-00 A^" 1 a-oc ^ p'(A n ) A"- 1 (A - A n ) ^ p'(A n ) 
We are done. 

(ii) From the definition of #'s for any k — 1, . . . , N — 1: 



?(A ) 

Using g(Afc) = p'{X k )p k we have 

It is easy to check that all p k > 0. There exists just one po such that p n = 1. This implies 
that #'s are coordinates and they are free. The Theorem is proved. 

The fact that the isospectral set S(Lq) is a diffeomorphic to R 1 ^^ 1 was already noted 
by Moser, [plfl . Tomei, [^7j] then showed that S(Lq) can be compactified and it becomes a 
convex polyhedron. The symplectic interpretation of this result as a version of the Atiyah- 
Guillemin-Sternberg convexity theorem was given by Bloch, Flaschka and Ratiu Q. 

Now for a function w(X) G Rat' N with Aq = we consider another, exponential represen- 



tation 3.2: 



w(X) = -^exp S(A), 



and define another set of angles by the formula 

Xk — X 



6' k = lim 



I (A) - S(0) + log 



A A: 

This formula can be put in the form 

N-l 



+ m k = 1, . . . ,N — 1. (5.5) 



7s Ik 



0' k = J u k + J u k + ni. 



' ' A s OO — 
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Whence, 9' k correspond to the Abel sum in the Baker form. The regularization is necessary, 
because the term J^ k uj^ diverges logarithmicaly on singular curve Y. Fathermore, we have 
simple relation 

e k = 6' k + m(k-l) + \og J] -^-A (5.6) 

The angles 9 and 9' differ by the real quantity which depend on the curve only. Evidently 
the variables 9' are coordinates on S(Lq) and their range is R 1 ^" 1 . 



6 The Poisson Bracket on Weyl Functions. 

Following [15| we consider functions w(X) with the properties i. analytic in the half-planes 



> and ^sz < 0. ii. w(z) = w(z), if ^sz ^ Hi. ^sw(z) > 0, if $sz > 0. All 
such function are called i?-functions. They play central role in the theory of resolvent of 
selfadjoint operators. The Weyl function of a Jacobi matrix is an i2-f unction. 

Atiyah and Hitchin introduced a Poisson structure on the space of rational functions. 
In the recent paper [[]] Faybusovich and Gehtman wrote Atiyah-Hitchin and higher Poisson 
structures on rational functions in compact invariant form. They defined the Atiyah-Hitchin 
bracket by the formula 

r / \ (W(A) — W(ll)) 2 

{w A ),w M } = \£1L. 6.1 

Here we discuss some of its' remarkable properties, [ 28| . 

We think about if (A) as an element of some commutative complex algebra which depends 
holomorphicaly on the parameter A. Evidently, |6.1| is skew-symmetric with respect to A and 
jj,. It is natural require linearity of the bracket 

{aw{\) + bw(X),w(u)} = a{w{\), w{v)} + b{w(X), w(v)}, (6.2) 

where a and b are constants. The symbol w(\) for A inside the contour C is given by the 
Cauchy formula 

1 ^ 2m J 



-d(. 



Whence due to |6.2| the values of the bracket in different points are related 



c 



{w(X),w( f i)} = — J 1 d C 



c 



It can be verified that the bracket BTT] satisfies this compatibility condition. 



Also it is naturally to require for the bracket the Leibnitz rule 

{w(\)w(fi),w(v)} = w(\){w([i), w(v)} + w(fj,){w(\),w(v)}. (6.3) 



preliminary draft 



15 



February 7, 2008 14:08 



It can be verified in a long but simple calculation that |6.2| - |6.3| imply the Jacobi identity 
{w(X), {w(/i),w(u)}} + {w(p), {w(u),w(X)}} + {w(y), {w{\),w{fi)}} = 0. 



The particular useful to us is an invariance of |6J] under the group of linear-fractional 
t r ansf or mat ions 

, aw + b ,„ 
w — > w = -, (0.4) 

cw + d v ' 

where a, b, c, d are constants. This property will be used in the construction of the second 
system of canonical coordinates. 

In our study of finite Jacobi matrices we need a small subclass Rat^ C R. These are the 
functions given by the formula |5.1| . All such functions have asymptotic expansion at infinity 

w(X) = - J ---.... 

We consider submanifold Rat' N with sq = 1, or equivalently, YlPk = !• We will demonstrate 
that the Dirac restriction of the bracket |0] to this submanifold takes the form 

{w(\),w(»)Y = HA) - wW) ( W(A ]~^ M - w(X)w(^ ■ (6.5) 

The linear Poisson structure on the phase space C is defined by the formulae 

{c k ,v k } = -Cfc/2, {c k ,v k+1 } = Cfc/2, (6.6) 

and all other brackets vanish. Whence, the linear bracket FTB corresponds to the restriction 



6.5 of the AH structure on the submanifold Rat' N . 

The formula |T| can be used to define the Poisson structure |BT6|. For example, substituting 



2.20 into |6.5|, after simple algebra we obtain 



X 2 jj? \fi X J A 2 /i 2 \fj, X 

From this one can read the first identity: {c ,v } = — c /2. 



A construction of canonical coordinates for the bracket |67T| or |6.5| will be given in terms 
of various representations for Rat^ and Rat' N . 



Canonical coordinates on RcLtjy. 
tion. 



The Dirac reduc- 



We start with the construction of the first system of canonical coordinate on Rat^ for the 
bracket IBTl]. The next Theorem shows that the parameters 



Xq, . . . , Ajv-i; Pqj ■ ■ ■ , Pn-i 



in the formula 5.1 are "almost" canonically paired. 
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Theorem 2 The bracket 6A in A — p coordinates has the form 

{Pk,Pn} = T r (l-5 fc ), 

{Pfc, A n } = Pfc^fc, 

{Afc, A n } = 0. 

Proof. We represent p's and A's as contour integrals 

1 



(7.1) 

(7.2) 
(7.3) 



Pk 



Pk A* 



2m 



2m 



o k 



In both formulae the contour Ok surrounding X k is traversed counterclockwise. Therefore, 
for k n 



{Pk,Pn} 



(w(C) -w(r])) 2 



(2m) 



C-v 



Wi 1 *' 



-dC, dr] 
1 



(2mfJ Jo n C-17 



d( dr] 



+ 



2m 



dr]w 2 (r]) 



On 



2m 



d( 



o k 



The two terms with square brackets vanish and for |7.4| we obtain 



2m 



d(w(() 



2m J Q-T) 



For ( in the exterior of the contour O p we have 



2m 



dr] 



o v 



w(rj) 



Pp 



(7.4) 



(7.5) 



Applying this formula twice to 7.5 we obtain |7.1|. If k — n then similar arguments show that 



the bracket 7.1 vanishes 



To prove |7.2j we compute for k ^ n 

{AfcPfc, p n } 
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From another side, 

{AfcPfc, Pn} = X k {p k , Pn} + Pk{Xk, Pn} ■ 

This together with imply that the bracket {X k , p n } vanish. For k = n we have 
{X k p k , Pk} = -Pk and {Pk, = Pk- The formula [T72| is proved. 

To prove |7.3| for k 7^ n we compute 

r , , -, 2X k p k X n p n 
lAfcPfc, ^npn) — — : : — • 

From another side, 

{Xkpk, X n p n } = PkPn{Xk, A n } + p k X n {X k , Pn} + PnXk{pki X n } + X k X n {p k , p n }. 

This together with formulae [Q] - [r3| imply that the bracket {Afc, A„} vanishes. For k = n 
arguments are the same. The proof is finished. 

Therefore, for the bracket |6.1| on Rat^ canonical coordinates are associated with poles 
of w{\): 

A , • • • , Aat__i, g(A ), . . • , q{\N-i)- 

Indeed, from the identity 

g(A) = \ - Pn 

p(X) ^Xn-X 



and condition p(X n ) = we obtain q(Xk) = p'{X k )p k . Furthermore, using fT2| - |P 

{g(A fc ), A„} = {p'(X k )p k , X n } = p'(X k )p k 5l = q(X k )S k . 
All other brackets vanish: 

{q(X k ),q(X n )} = {A fc ,A„} = 0. 

In this coordinate form the Poisson structure on RatN was introduced in ||. These identities 
imply 

{q(X),q(p)} = {p(X),p(p)} = 

and 

q(X)p(fi) - q(p)p(X) 



X — p 



The last expression is called a Bezoutian, see [TBI]. This form of the bracket easily leads to 



6.1 



The second set of canonical coordinates on Rat^ is associated with zeros of w(X): 

7i,...7at-i; p(7i),---,K7at-i); qo,Po- 



To prove this, we introduce the new "dual" function w'(X) as an inverse of the function |572 

1 p(A) 



w'(X) 



w(X) q(X) 
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Due to O the bracket for the dual function is given by the formula xT. The new meromor- 
phic function maps the upper half-plane into itself and has the expansion 



w 



. N-l , 

'(A) = - + c+ V^-, 



where, 



PoQo - qi ft , E7 S , ^ n 

c = s = 1 , p s > 0. 



4 



go go 



Theorem 3 The following identities hold 



{pk'Pn} 


2 P'kPn ^ 
7n - 7fc 




= M 


{7fc,7s} 


= o, 


{go,Pfc} 


= {go,7s} = o, 


{p'k,Po} 


= Pi 




= o, 


{po,go} 


= go- 



Proof. The identities |7.6| - [7.8| using integral representation 

1 



Pk 



Pk Ik 



2m 
1 



2m 



C«/(C) 



can be proved exactly the same way as identities |7.1| - [7.3| of Theorem 2. 



Let us compute the first bracket [FJ5 

{qo,p'k} ■ 



A 



i 



w'(()d(} 



A->oo io'(A) ' 2m' j 0k 

lim A f (w'(X)-w'(Q) 2 
x^oo 2mw'{\) 2 J 0k A - C 



d(. 



(7.6) 

(7.7) 
(7.8) 
(7.9) 
(7.10) 
(7.11) 
(7.12) 



Expanding the square and computing each term separately we see that the bracket vanishes. 
The proof of the second identity [T9| is exactly the same. 

To prove |7.10| we note 

{Po» Pk) _ P± 

go go' 



{c, Pk) 
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From another hand 



{c,Pfe} = {lim w'(X) - K-^-. 



w'(C)d(} 



o k 



lim 

a->oo 2ni 



o k 



X(A) - w'{Qf 
A-C 



Comparing it with the previous formula we obtain 77T0 . The proof of the formula [7.11| is 
similar. 



1. 



To prove the last formula [7.121 we compute 

{POlQo} r -i r ,- //\n A 

= {c,qo\ = { hm w (A) ,g | 

= hm{w (A), hm } = hm hm — — 

A^oo V-kx> W' {[!) A-+ooft->oo w'{jJ,) z A — [I 

This implies the result. The proof is finished. 

From [771] using the formula p(7 s ) = —p ! s q\ls) we obtain 

From [7g, [777| and [73 

{p(7n),P(7fc)} = 0. 

From [7j9 

{p(7n),go} = o. 

From |77ig and [TTTg 

{p(7n),Po} = 0. 

These identities together with identities of Theorem 3 provide a proof of our statement. 

This coordinate system is usefull in construction of the Dirac restriction, M, of the AH 
bracket |67J on the submanifold M C -Rat at determined by the conditions 

$i = Po = ci, $ 2 = log g = c 2 

where c\ and c 2 are some real constants. 

Consider a more general problem. The submanifold M of dimension 2N—m is determined 
by the conditions 



$i = Ci, $ 2 = c 2 , 



where $'s are some functions on the phase space and c's are real constants. The bracket 
{•, •} is modified 



{F u F 2 }' = {F u F 2 } + *k{Fi, ®k} 



fe=i 
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whith cr's chosen such that 



$l = {$ k ,F 2 } + J2<7s{®k,$s} = 0. 



8=1 



for all k — 1, . . . , m. Geometricaly, this condition means that the vector fields produced in 
the bracket {•,•}' are tangent to M. If the matrix ||{$fc,$ s }|| has an inverse ||Cfc a ||, then 
the last system can be solved for cr's and 

{F 1 , F 2 }' = {F x , F 2 } +J2i F ^ C *> {^2, $ s }- 

k,s=l 

Implementing this procedure for our choice of functionals $i and $2 we obtain 

o x = {log g , ^2}, 0-2 = -{^0,^2} 

and 

{F 1 ,F 2 }' = {Fi,F 2 } + {\ogq ,F 2 }{F uPo } - {p , F 2 }{F 1; logg }. 
One can easily verify 

{?(A),?(A*)} , = fr(A),p(A*)} / = 0. 
Using {po,<?(A)} = <?(A) and {p(p),qo} = q(p) we obtain 

{ q (\), P (»)Y = q{X)M " qifl)p{X) + g(/i)g(A) 



Finaly, we have 



{w(X),w( t x)}'^(w(X)-w( f i)) 



A - 11 q 

w(X) — w(\)w(fj,) 



A — p 



This becomes |6.5| for a particular choice $2 = logc/o = 0. 



8 The canonical coordinates on Rcit'jy. 



Now we turn to the submanifold Rat' N with q = ^2 Pk — 1 and the Poisson bracket |T5 



Here the situation is a little more subtle. In all formulae we omit the prime near the bracket 
{•,•}'• 



Theorem 4 The bracket \6.3{ in A — p coordinates has the form 

2pk Pn I sr^ Ps Ps 



{pk,pn} 



An. — Al 



ZpkPn 



X s — Afc A s — A n 



{PkAn} = ~PkPn + Pk^k, 
{Afe, A n } = 0. 
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Proof is similar to the proof of Theorem 2 and therefore is omitted. 
The Theorem implies, 

{g(A fc ), A„} = -q(X k ) Pn + q(X k K. (8.4) 
Thus, we have the first system of canonical coordinates, so-called, the action-angle variables 

Ai, . . . , Ajv-i; Oi, ■ ■ ■ , On-i (8-5) 
where 0's defined by [5.3| are real and canonically paired with A's. Indeed, |S.4| implies 



The formulae 8J] and KM produce 



{logg(A fc ), logg(A„)} = J2 f^r " £ f^T- ( 8 - 6 ) 



This identity implies the commutativity of the angles: {9k, 9 n } = 0. 

Using the Theorem it can be checked easily that po = — A& is a Casimir of the bracket 



Evidently, for the restricted bracket the canonical relations established in Theorem 3 
survive. We have the second set of canonical variables on Rat' N 



7i,...,7jv-i,7ri,...,7rAr-i; 7r fe = log(-l) 7V+ V(7fc)- (8.7) 



These divisor-quasimomentum coordinates were introduced in ||T8|. The denominator p(X) 
(~ l ) N E[(An - A) satisfies 

(-l) N+k p( lk ) > 0, k = l,...,N-l. 

Whence 7r's are real and canonicaly paired with j's: 

{TT n ,Jk} = Sl 

All other brackets vanish. 



In the rest of this section we show that the variables 5.5 associated with representation 



of w(X) G Rat' N in the exponential form |372 



w (\) = -\e^ x) 
A 

can be moved by corresponding A's: 

{9' k ,\ n } = S k n , fc,n = l,...,iV-l. 

Though it follows from the previous discussion of the action-angle variables and the formula 
|5.6| we will give an independent proof of this fact. It is important to notice that we can not 
expect commutativity of the variables 9' n . 
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The multivalued function H(A) has the form 5(A) = ^^T^ 1 log(7 s — A) — log(A s — A) 
and defined up to an integer multiple of 2iri. The bracket { • , 5(A)} is single valued since 
additive constant vanishes. The Poisson bracket |5.5| in terms of the function 5(A) has the 
form 



{2(A), S(/i)} 



4sinh^(S(A) -E(ji) - log A + log/i/2) l s 



A — n 



+ _ e =(A) _ ± e HM 
A 



or 



{S(A),E(/i)} 



1 



- w(X) + w{jj) 



w(X)w(p) A — ji 

which is more convinient for calculations. The pole Xk can be represented as a contour 
integral 



o k 



Or, 



where X' k is an arbitrary fixed point on the contour Ok surrounding A&. As a simple example 
we prove commutativity of A's 



{Afc ' A4= (2^/ J{ E (0,Z(v)}d(d V 

O k On 



From the formula R7SI one can easily see that the double integral vanishes. 
Now for the angles 9' n we have 

{6' k ,X n } = lim [{2(A), A n }- {2(0), A n } + {log A fc -A/A fc ,A n }] 

The last term vanishes. The first two terms are more complicated 

1 



{2(A), A n } 



2iti 



O n 



{~(X),~(()}d( 

w(X) d( 
2m J w{()\-( 



0„ 



d( 



2m J A-C 

O n 



+ 



2ni J w(X) A-C 



O n 



1 

2tH 



(8.9) 
(8.10) 
(8.11) 
(8.12) 



O n 



If A — > Afc, k n, or A — > Ao = 0, the terms |8.9| - |3.11| vanish. The term |8.12| is equal to 
— p n . Therefore, 



lim {5(A), A„} = -p r 

A— >Afe 
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This implies {9' k , X n } = 0, k ^ n. 



Furthermore, if A 



X n , then the term 8.S becomes —1, the term 8.10 becomes 2 and 



S.ll vanishes. The term 8.12 is — p n . Therefore, 



lim {E(A), A n } = I- p n 

A— »A n 



Thus {9' n ,X n }=l. 



9 Tangent and Transversal Flows. 

Using the poles of w(X) we define Hamiltonians Hj — ~ J] A^, j = 1, . . . , N. The flows 



produced by them in the bracket |T5] are tangent to the isospectral manifold: {A*,, Hj} = 0. 



Due to Theorem 4 the standard, [^T[j, Toda flows have the form 

p' k = {p k , Hj} = (Af 1 - K~ X Pn) Pk, k = 0, . . . , N - 1. 

Toda flows commute {Hj, H^} = and linearized in the variables |B.5| : 

61 = {6 k , Hj} = (XI 1 - XI 1 ), k = 1, . . . , N - 1. 

Similarly, from zeros of w(X) we define another set of Hamiltonias 2} = j^'yl, j = 
1, . . . , N — 1. By Theorem 3 the flows produced by these Hamiltonians do not affect 7's. 
Therefore, we call these commuting flows transversal. They are linearized in the variables 

nl = {n k ,T. j } = 1 l 1 , k = l,...,N-l. 



This is an example of the situation similar to the one considered in physics, |TT| , |J 



Given two systems of canonical coordinates and two families of commuting Hamiltonians. 
Each family depends only on the half of coordinates of the corresponding canonical system. 
Hamiltonians of both families produce coordinate system for the Poisson manifold. 

It is routine exercise to derive the equations of motion for 7's under H flows and for 
A's under T flows. Then the inverse spectral problem can be solved using trace formulae of 
section 2. We do not dwell on this. 
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